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摘要 : 围绕 着 恒星 级 黑洞 的 中 微 子 主导 吸 积 盘 可 以 通过 盘 上 发 出 的 中 微 子 淫 灭 为 件 玛 暴 
提供 能 量 。 对 于 黑洞 超 吸 积 系统 ， 吸 积 可 能 引起 黑洞 特征 的 极 大 演化 ， 而 这 会 进一步 引 
起 中 微 子 光度 的 演化 。 考 虑 不 一 样 的 平均 吸 积 率 和 初始 黑洞 参数 ， 通 过 吸 积 系统 的 演化 
分 析 中 微 子 潭 灭 光度 和 总 的 中 微 子 削 灭 能 量 随时 间 的 变化 。 同 时 计算 了 短 暴 GRB 
090510 的 中 微 子 淫 灭 能 量 并 且 与 理论 预测 的 结果 比 对 ， 发 现 NDAF 的 中 微 子 淫 灭 能 量 


的 理论 预测 值 远 高 于 观测 值 ， 意 味 着 这 种 模型 可 能 可 以 提供 GRB 090510 爆发 所 需要 的 


能 量 。 
关键 词 : PX. USE. 黑洞 ， 伽 玛 暴 
中 图 分 类 号 : P172.3 文献 标识 码 : A 文章 编号 : 1672-7673 (2018) 


伽 玛 射线 暴 (Gamma-Ray Bursts, GRBs) 是 来 自 宇宙 空间 的 伽 玛 射线 在 短 时 间 内 
忽然 增强 的 高 能 爆发 现象 。 它 的 持续 时 间 一 般 在 0.1s 到 1000s 之 间 ， 通 常用 Too (Tooti 
光子 数 累 积 计 数 从 5% 到 95%) 作为 典型 时 标 ， 把 伽 玛 射 线 暴 分 为 长 暴 (T%>2s) 和 短 暴 
(Teo<2s) 两 种 类 型 .两 类 暴 可 能 产生 于 不 同 的 辐射 区 域 , 它们 的 暴 源 可 能 有 本 质 的 区 别 凸 。 
目前 普遍 认为 长 暴 起 源 于 大 质量 恒星 的 引力 塌 缩 ， 短 暴 来 源 于 双 中 子 星 或 者 中 子 星 与 黑 
洞 双星 系统 的 并 合 。 引 力 塌 缩 或 者 双星 并 合 都 可 能 形成 一 个 以 恒星 级 黑洞 为 中 心 ， 周围 
是 一 个 吸 积 盘 的 超 吸 积 系统 中。 在 吸 积 系统 中 ， 通 常 可 以 认为 是 中 微 子 潭 灭 过 程 为 伽 玛 
暴 提供 能 量 ， 而 一 个 围绕 着 恒星 级 黑洞 的 具有 非常 高 的 吸 积 率 的 中 微 子 主导 吸 积 才 

(Neutrino Dominated Accretion Flow, NDAF) 被 认为 是 伽 玛 暴 中 心 能 源 机 制 的 重要 候 
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选 者 SH。 中 微 子 主导 吸 积 盘 的 内 区 吸 积 物质 具有 非常 高 的 密度 和 温度 ， 质 子 和 光子 无 
法 从 盘 上 逃逸 。 相反 , 中 微 子 可 以 带 着 粘 灌 产生 的 热量 和 大 量 的 引力 束缚 能 从 盘 上 辐射 ， 
而 它们 在 盘 上 的 漂 灭 可 能 产生 初始 火球 。 

文 [5] 推 导 了 通过 中 微 子 提供 能 量 的 喷 流 的 最 大 可 能 能 量 , 该 能 量 是 暴 持续 时 间 的 函 
数 ,结果 证 明了 这 个 模型 可 以 为 持续 时 间 小 于 100s 的 伽 玛 暴 提供 能 量 ， 而 晚期 的 X 射线 
耀 发 和 超 长 伽 玛 暴 很 难 用 这 个 模型 解释 。 但 是 ， 他 们 忽略 了 来 自 于 中 心机 制 的 活动 ， 以 
及 黑洞 的 演化 ， 特 别 是 黑洞 转速 的 演化 。 在 中 微 子 主导 吸 积 盘 中 ， 中 微 子 辐 射 和 中 微 子 
潭 灭 光度 与 吸 积 率 、 黑 洞 质量 以 及 黑洞 转速 紧密 相关 负 。 如 果 在 伽 玛 暴 的 中 心 ， 真 的 存 
在 一 个 围绕 着 中 微 子 主导 吸 积 盘 的 恒星 级 黑洞 ， 那 么 在 这 种 黑洞 超 吸 积 系统 中 ， 吸 积 硬 
能 引起 黑洞 的 剧烈 演化 ， 黑 洞 参数 不 可 避免 地 发 生 巨大 变化 ， 而 这 将 进一步 影响 中 微 子 
潭 灭 光度 以 及 总 的 中 微 子 能 量 中。 首先 讨论 不 同 的 初始 参数 下 转速 的 演化 ， 肯 定 了 黑洞 
的 演化 引起 黑洞 参数 的 巨大 变化 。 接 着 讨论 不 同 参数 下 中 微 子 光度 和 总 的 中 微 子 能 量 的 
演化 ， 分 析 黑 洞 参数 对 中 微 子 光度 和 能 量 的 影响 。 最 后 根据 短 暴 GRB 090510 的 观测 数 
据 计 算 GRB 090510 的 中 微 子 能 量 ， 并 与 理论 预测 值 进行 比较 ， 分 析 这 种 机 制 解 释 短 暴 
GRB 090510 爆发 的 可 行 性 。 


1、 物 理 模 型 分 析 
理论 上 ， 如 果 在 伽 玛 暴 的 中 心 存在 一 个 超 吸 积 系统 ， 那 么 黑洞 的 特征 参数 将 随 着 时 
间 发 生 很 大 的 变化 。 不 考虑 BZ 机 制 ， 根 据 能 量 和 角 动 量 守恒 公式 ， 可 以 得 到 一 个 Kerr 


黑洞 的 演化 方程 中: 
QM py _ Me ， (D 
dt 
Wan -MI ， (2) 
dt 


其 中 Mi > Jay 41M 分 别 是 黑洞 的 质量 、 角 动量 和 平均 吸 积 率 。 下 面 用 无 量 纲 化 的 参数 


n= 1 2 M ^ Al E 
my, =M gy I MUR = MIM oS) page, Em 和 im 对 应 于 最 大 稳定 半径 Tw 处 的 能 量 和 角 
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Reb, 4 mola /GM 为 黑洞 的 无 量 纲 转速 ，Yw = 342. -VB-Z)B+Z, +22) 为 无 


纲 化 盘 的 最 大 稳定 轨道 半径 。 其 中 ，21=1+0-a) I0) «07 a)^LZ, = Ba? +Z. 
根据 (1) - (4) 式 可 以 得 到 黑洞 转速 随时 间 的 演化 : 
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如 果 黑 洞 初始 的 质量 四 、 转 速 w 和 吸 积 率 六 给 定 ， 那 么 根据 方程 可 以 得 到 任意 时 
刻 黑洞 的 特征 参数 。 而 吸 积 率 % 的 给 定 ， 则 需要 考虑 吸 积 紊 和 盘 质 量 之 间 的 关系 ， 盘 质 


量 可 以 根据 下 式 估算 ; 
M, =M 22 
I+z , (6) 


FLA Too 和 z 分 别 是 暴 的 持续 时 间 和 红 移 。 为 了 简单 方便 ,用 无 量 纲 化 的 参数 
Maa =M ask /Me 计算 。 


在 中 微 子 主导 吸 积 盘 模型 中 ， 由 于 中 微 子 潭 灭 光度 L,; 与 黑洞 质量 MM, 、 黑 洞 转速 


a 和 吸 积 率 区 有 密切 的 关系 ， 所 以 黑洞 质量 和 角 动 量 的 剧烈 演化 也 会 引起 中 微 子 光 度 的 
演化 。 中 微 子 光度 的 理论 分 析 表 达 式 名 
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供 ， 因 此 ， 中 微 子 漂 灭 总 能 量 E,; 等 于 根据 瞬时 辐射 阶段 的 观测 数据 得 到 的 各 向 同性 辐 


射 能 量 人 7 和 长 时 标 余辉 能 量 外 流 的 各 向 同性 动能 By 之 和 D3l。 
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其 中 ,7 为 转化 效率 ;2 为 喷 流 的 张 角 。 选 取 刀 = 0.3 ， 可 以 从 观测 数据 计算 出 中 微 子 总 
的 潭 灭 能 量 。 


2 结果 和 分 析 
根据 上 面 的 方程 ， 可 以 通过 合理 地 给 定 黑洞 的 初始 质量 m PR a, FUR m , 
得 到 中 微 子 潭 灭 光度 和 总 的 中 微 子 潭 灭 能 量 随时 间 的 变化 。 
参数 选择 的 合理 性 是 一 个 至 关 重要 的 问题 ， 而 盘 质量 是 反应 这 种 模型 合理 性 的 一 个 
重要 因素 。 短 暴 来 源 于 致密 双星 的 并 合 ,因此 ,黑洞 的 质量 通常 小 于 双星 的 总 质量 , 即 ; 
< 4M。， 双 中 子 星 或 者 黑洞 与 中 子 星 双星 并 合 的 模拟 结果 表明 , 盘 质 量 大 约 在 0.2M。 到 


pui 


0.5M o ZEPS 可 以 合理 地 假设 黑洞 初始 质量 m, = 2.3413, AE Ema, = 0.2 和 0.5, 选取 


黑洞 转速 a = 0.0.5 和 0.9。 根 据 CIO AE G) 式 得 到 转速 a 随时 间 变 化 的 关系 曲线 ， 


图 1( 左 ) 黑洞 转速 a 的 演化 。 实 线 和 虚线 分 别 表示 初始 黑洞 质量 为 2.3 和 3. 


图 2〈 右 ) 中 微 子 漂 灭 光度 工 :的 演化 。 所 有 线 的 初始 黑洞 质量 取 2.3. 实 线 和 虚线 分 别 
表示 吸 积 率 为 0.2， 初 始 转速 为 0.5 和 0.9. 虚 点 线 和 点 线 分 别 表 示 吸 积 率 为 0.5， 初始 转 
速 为 0.5 和 0.9. 


Fig.1(Left) Evolution of the spin of BH a. Solid and dashed lines correspond to the initial 
dimensionless value of BH mass 2.3 and 3, respectively. 


Fig. 2 (Right) Evolution of the mean neutrino annihilation luminosity L,.. All lines 


correspond to the initial dimensionless value of BH mass 2.3. The solid and dashed lines 
correspond to the dimensionless mean accretion rate 0.2 and the initial spin of 0.5 and 0.9, 
respectively. The dashed dotted and dotted lines correspond to the dimensionless mean 
accretion rate 0.5 and the initial spin of 0.5 and 0.9, respectively. 
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如 图 1。 图 1 中 ， 初 始 转速 分 别 取 0 和 0.9， 实 线 和 虚线 分 别 表示 初始 黑洞 质量 为 2.3 和 
3. 从 图 中 可 以 清楚 地 看 到 黑洞 转速 随时 间 的 演化 有 比较 明显 的 变化 。 而 黑洞 的 转速 对 
吸 积 盘 质 量 的 影响 尤为 重要 ， 因 此 在 计算 中 必须 考虑 转速 的 演化 。 从 图 中 也 可 以 看 出 黑 
洞 初始 质量 对 转速 的 演化 影响 比较 小 ， 所 以 后 面 的 计算 中 ， 统 一 选取 黑洞 的 初始 质量 为 
2.3。 接 着 根据 〈7) 式 计算 中 微 子 淹 灭 光度 的 演化 ， 如 图 2， 图 中 实 线 和 虚线 分 别 表 示 
吸 积 率 为 0.2， 初 始 转速 为 0.5 和 0.9 的 情况 ， 虚 点 线 和 点 线 分别 表 示 吸 积 率 为 0.5， 初 
始 转速 为 0.5 和 0.9 的 情况 。 从 中 可 以 看 到 中 微 子 光 度 随 着 时 间 的 演化 变化 很 大 ， 它 的 
最 大 值 可 以 达到 10°’erg s! 

中 微 子 潭 灭 光 度 的 计算 与 观测 数据 不 具有 直接 的 相关 性 。 为 了 进一步 检验 中 微 子 漂 
灭 模型 的 合理 性 ， 选 择 短 暴 GRB 090510 的 观测 数据 四， 如 表 1， 计 算 了 短 暴 090510 的 
中 微 子 潭 灭 总 能 量 E,; ， 并 且 与 理论 计算 分 析 线 比较 ， 如 图 3。 根据 (6) 式 ， 选 取 GRB 


090510 的 Too Al z, img =0.2 时 ， 吸 积 率 为 1.27。3 条 线 对 应 的 参数 分 别 是 虚线 : 吸 
积 率 1.27， 初 始 转速 0.9， 实 线 : 吸 积 率 1.27， 初 始 转速 0.5;， 虚 点 线 : 吸 积 率 0.2， 初 
始 转速 0.5。 星 星 表示 根据 GRB 090510 的 观测 数据 由 (9〉 式 计算 出 来 的 的 E,; ， 发 现 
GRB 090510 的 中 微 子 潭 灭 总 能 量 的 观测 值 远 小 于 理论 预测 值 ， 而 且 初 始 转 速 或 盘 质 量 
取 更 大 值 时 ， 吸 积 率 也 将 更 大 ， 这 时 ES 的 理论 预测 值 更 大 。 


表 1 GRB 090510 的 观测 数据 
Table ] Observational data of GRB 090510. 


GRB Too Z Fy FX(1 hours) Ey,iso Ex,iso 0 j E us 
(s) (10"erg cm?) (10 erg cm? s) (10“ergs) — (10lergs) (rad) 
(10“ergs) 
090510 0.3 0.903 3.4 5.04 0.73 MMS 0.017 0.028 
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t 
图 3: 中 微 子 总 能 量 理论 预测 值 和 观测 值 的 比较 。3 条 线 黑洞 初始 质量 都 为 2.3， 虚 点 线 
表示 吸 积 率 为 0.2, 初始 转速 为 0.5; 实 线 和 虚线 都 是 吸 积 率 为 1.27, 初始 转速 分 别 为 0.5 
和 0.9， 线 终止 的 位 置 为 转速 大 于 0.998 的 位 置 ， 星 星 表示 GRB 090510 的 观测 值 。 


Fig.3 Comparison of the theoretical predictions and observational date of the total 
annihilation energy. Three lines all correspond to the initial dimensionless value of BH mass 
2.3。 The dash dot line correspond to the dimensionless mean accretion rate 0.2 and the initial 
spin of 0.5; the solid and dashed lines correspond to the dimensionless mean accretion rate 
1.27 and the initial spin of 0.5 and 0.9, respectively. The lines truncate when the BH spin has 
evolved to 0.998. The star denote the observational result of GRB 090510. 

3 结论 和 讨论 

短 暴 最 有 可 能 起 源 于 双 致 密 星 并 合 ， 因 为 这 种 情形 产生 的 吸 积 盘 较 小 ， 持 续 时 间 可 
以 和 暴 的 持续 时 标 相 当 , 而 围绕 着 恒星 级 黑洞 的 中 微 子 主导 吸 积 盘 可 以 通过 盘 上 发 出 的 
中 微 子 潭 灭 为 伽 玛 暴 提供 能 量 。 对 于 这 样 一 个 黑洞 超 吸 积 系统 ， 盘 质量 是 反应 这 种 模型 
合理 性 的 一 个 重要 因素 ， 根 据 短 暴 的 数据 ， 在 火球 和 中 微 子 潭 灭 模型 的 框架 下 ， 用 黑洞 
的 不 同 参数 估算 了 吸 积 盘 的 质量 ， 发 现 短 暴 吸 积 盘 的 质量 主要 依赖 于 输出 能 量 、 喷 流 的 
张 角 和 黑洞 的 参数 。 而 黑洞 参数 中 ， 黑 洞 的 转速 和 吸 积 率 对 盘 质 量 的 影响 尤为 重要 中。 

选取 不 同 的 初始 黑洞 转速 和 平均 吸 积 率 ， 可 以 明显 地 看 到 转速 随时 间 的 变化 很 大 ， 
即 : 吸 积 可 能 引起 黑洞 特征 的 极 大 演化 ， 而 这 会 进一步 引起 中 微 子 光度 的 演变 。 进 一 步 
讨论 了 不 同 吸 积 率 ,给 定 不 同 的 黑洞 初始 参数 的 情况 下 ,平均 中 微 子 潭 灭 光度 的 含 时 演 
化 以 及 总 的 中 微 子 潭 灭 光 度 的 演化 。 为 进一步 讨论 这 种 模型 解释 伽 玛 暴 瞬 时 辐射 和 余辉 
阶段 的 可 能 性 ， 选 取 了 短 暴 GRB 090510 的 观测 数据 ， 通 过 观测 数据 计算 了 伽 玛 暴 观 测 
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的 中 微 子 潭 灭 总 光度 ， 并 且 和 上 面 的 理论 预测 结果 相 比 较 发 现 ，GRB 090510 的 中 微 子 
潭 灭 总 光度 在 理论 预测 线 下 方 ， 远 小 于 理论 预测 值 ， 这 可 能 意味 着 这 个 理论 模型 可 以 提 
供 GRB 090510 爆发 所 需要 的 能 量 。 
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The neutrino annihilation model of NDAF explains the explosion 
of GRB 090510 
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Abstract A neutrino-dominated accretion disk around a stellar-mass black hole can power a gamma-ray 
burst via annihilation of neutrinos launched from the disk. For the black hole hyperaccretion system, 
accretion should trigger the violent evolution of the black hole’s characteristics, which further leads to the 


evolution of the neutrino annihilation luminosity. In this paper, we analyze the mean time-dependent 
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neutrino annihilation luminosity and the total neutrino annihilation energy via the evolution of the accretion 
system by setting the different mean accretion rates and initial black hole parameter. Furthermore, we 
calculate the total neutrino annihilation energy of the short GRB 090510 by observational data, and 
compare with the results of theoretical analysis. It is found that the theoretical valre of the annihilation 
energy of NDAF neutrinos is much higher than that of the observed values, which means that the 


mechanism may can explain the explosion of the short GRB 090510. 
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